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Abstract 

We demonstrate that an n'^ /i/n'^ junction is the most suitable candidate for electronic control 
and readout of qubit states in quantum computing systems based on shallow impurities. The 
signature of this system is that the n''"— regions serve as metallic electrodes separated form the 
i— region by a self-induced barrier (internal workfunction) . The n+/i/n+ system mimics the prop- 
erties of a metal-vacuum-metal junction with the qubit (impurity atom) placed in a "vacuum" 
i-region between two "metallic" n"*" electrodes. We will show that the self-induced barrier exists in 
a sufficiently wide range of the concentration of dopants in the n"^-semiconductor (e.g. up to 10^^ 
cm~^ for Si) and its height can be controlled by tuning the doping level. A shallow donor placed 
in a vacuum i-region will be populated with one electron in equilibrium. In the case of Li donor 
in Si the n"'"-electrodes will be used for a precision placement of the Li atom during the growth 
process; for voltage control and manipulation of the qubit states; and for a qubit readout by means 
of the optically stimulated resonant tunnelling. Another important feature of our system is that 
the qubit states (first two lowest energy levels of Li in Si) are separated by an energy gap from a 
continuum of the many-body states of the controlling electrodes. 
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Introduction. Various systems have been proposed during last decade as potential candi- 
dates for solid state quantum computing (QC). All these systems are based on control of 



either localized spins or localized charged states in semiconductors. In particu 



proposed to encode qubits on: nuclear spins of '^^P donors or Si isotopes l3|| m bi ma- 

iiiy,Qy,y,Q; 



ar, it was 
in Si ma- 



trix; spins of electrons localized in semiconductor quantum dots (QDs) 
and on charged states of shallow impurities in Si 111 . For all these cases the distance 
d between neighboring qubits ranges from 5nm ^ J^, Ij to lOOnm 11, Q] and electronic 
control must be used for a readout and for one- and two-qubit gates. It is clear, that the 
distance a between an individual qubit and its controlling electrode has to be smaller than 
the interqubit separation d to ensure a proper space resolution of the qubits. Thus the 
distance a must be in the range a < 5nm (nuclear spins [1,^0]) and a ~ 20nm (shallow 
impurities in Si |ll[ [l^ ) . 

It seems almost obvious that one can use metal electrodes for control and readout of the 
solid state qubits. As it has been suggested in most of the solid state QC proposals 



ll| , the thin-film metal electrodes have to be deposited directly on a Si surface or on a thin 
insulating Si02 layer grown on this surface. A careful analysis of this idea reveals that a 
single shallow donor in an undoped semiconductor placed in a proximity (a ~ 5-20nm) of an 
interface with the metal will be ionized. It simply means that the qubit cannot be formed. 
This conclusion is in agreement with the fact the localized electrons on shallow impurities 
placed in vicinity of metal electrodes have not been observed in QC experiments. Indeed, it 
is well known that a high potential barrier (a Schottky barrier) is always formed at metal- 
semiconductor (MS) or metal- insulator-semiconductor (MIS) interfaces The barrier 
height is determined by a large density of the interfacial states that pin the Fermi level Ep 
at energy A ~ 0.5-0.8 eV below the conduction band edge Ec{0) at the interface. The pinning 
occurs for Si and GaAs in contacts with almost all metals and for doping concentrations 
N < 10^^ cm~'^ the barrier width la is in the range Id ^ 100 nm However an extremely 
pure (i-)S layer with defect density < 10^° cm~'^ has to be used for the solid state QC. 
For these concentrations the barrier width 1^ is in the range of /im or more, and inside iS 
layer the conduction band edge Ec{z) is practically constant at the distances to the interface 
corresponding to the qubit location, 2; < a ~ 20 nm. Therefore the potential of the i-S layer 
between two control metal electrodes will be flat (see Fig. ^a)). The barrier height is in 
the range 0.5-0.8 eV and all shallow donors with ionization energy Id < 100 meV will be 



ionized. 

One can apply an electric field to push the impurity level below the Fermi energy of the 
left electrode (see Fig. Hl^a)). However the required electric field, E ~ A/ae > 3 x lO^V/cm 
(e is the elementary charge), is so strong that the impurity will be ionized due to the Zener 
tunnelling breakdown into the right contact isl ]. 

The potential barrier can be efficiently decreased if a very thin heavily doped n'^-S layer 
with donor concentration ~ 10^° cm~^ is formed between M and i-S layers The 
required conditions will be fulfilled if Ep < Ec{a) < Ep + Id in ^-semiconductor layer at a 
distance a (cf. Fig. ^a)). Unfortunately, these extremely stringent conditions cannot be 
satisfied simultaneously for a sufficiently large interface plane and a sequence of donor qubits 
located near it, because of the large fluctuations of the Schottky barrier A and the number 
of dono. ,„ ..e .a..tlQy Q. Sin,,,. ove,wM.„„ ob..,e. ..o .i. 

in MIS structures with thin insulating lay ers (Fig. ^b)). We note that the barrier height 
Ai at Si02-Si interface is about 4 eV |13] and therefore it is impossible to push the donor 
level below Ep with any realistic electric fields. 

The purpose of this paper is to demonstrate that an /i/n^ junction is the most suitable 
candidate for electronic control and readout of qubit states in quantum computing systems 
based on shallow impurities. The signature of this system is that the n"*"— regions serve as 
metallic electrodes separated form the z— region by a self-induced barrier (internal work- 
function). The jijn^ system mimics the properties of a metal- vacuum-metal junction 
with the qubit (impurity atom) placed in a "vacuum" i-region between two "metallic" n"*" 
electrodes. As we will see below, the self-induced barrier exists in a sufficiently wide range 
of the concentration of dopants in the n^-semiconductor (e.g. up to 10^° cm~^ for Si) and 
its height can be controlled by tuning the doping level. A shallow donor placed in a vacuum 
z-region will have its ground state energy below the Fermi level and will be populated with 
one electron in equilibrium. 



121 the same n"'"-electrodes 



Finally, we note that for a qubit based on Li donor in Si 
can be be used for a precision placement of the electro-migrating Li ions during the growth 
process; for voltage control and manipulation of the qubit states in Li at low temperatures; 
and for a qubit readout by means of the optically stimulated resonant tunnelling. 



LDA treatment of doped semiconductor "'"-i structure. Let us consider formation 




FIG. 1: Schematic band diagrams of Metal-Si-Metal (a), and Metal-Si02-Si-Metal junctions (b). 



of a barrier $s at a ri^ /i interface. We will assume that the donor concentration in rt^- 
region exceeds the threshold of the metal-insulator phase transition. If we replace spatial 
distribution of ionized donors with positive uniform charge and consider a single-valley 
parabolic-band n-type semiconductor (e.g . GaAs) the problem can be mapped onto the 
problem of a metal- vacuum interface jisl Ifol I20I . The barrier of an ri^ /i junction or, 
equivalently, the workfunction of a simple metal can be expressed in the form known as 



Budd-Vannimenus theorem 
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[oo)-0(O)-£. (1) 

The meaning of Eq. (PJ) is illustrated in Fig. |21 where a schematic band diagram of the 
n+/i-junction is shown. The energies E'^ and E^j are the conduction band edges of ri^ and 




FIG. 2: Band diagram of /i junction and schematic behavior of (/'(x) and Vcs{z) 



^-semiconductors respectively, V"cfT(-2) is the exact one-electron quasi-particle potential, and 
(f){z) is the exact electrostatic (Hartree) energy of the electron related to the electron density 
n{z) through the Poisson's equation: 



(2) 



where Q{y) is a step function. The quantities ft and e are the exact chemical potential 
many-body effects included [2^ and the exact total energy per one electron in the bulk 
respectively. All these quantities are well known within the local density approximation 
(LDA). In particular, 

5 4 \ r| Ts + 7.S J 

where we introduced the effective Rydberg: 

2xao 2m* Og 

X is the dielectric constant, m* is the effective mass, and is the effective Bohr radius: 



Ry* = 7^ = (4) 



«o = — r^- (5) 

me 

The parameter is related to the bulk electron concentration no in the region as follows: 

1/3 / Q \ 1/3 

(6) 



47^000/ \47rA^ao 
where N is the donor concentration and we assume that all the donors are ionized, i.e. 
no = N. 

Since e is known the problem of finding the barrier $s is reduced to calculation of the 
potential difference (pioo) — 0(0) caused by the interfacial dipole layer. The latter can be 
expressed as: 



0(oo)-0(O) = — '^^^^V^/" ^ ,\ (7) 



f^g{e)deJ^-mC)fdC 
where we introduced dimensionless energy e = [E — Vcsi—00)]/ Ep, and coordinate ( = 
kpz = (^tt'^NY^^z. Here Q is the total number of electrons per unit area in the half-space 
C < 0, and 

/ X 3 1 — e 

= J-^ (8) 

is a dimensionless effective density of the occupied states for electrons moving perpendicular 
to the interface. 

Eq. (|7j) is exact and follows from the properties of the 1-dimensional Poisson's equa- 
tion. The internal workfunction $s can be calculated exactly provided that we know the 
eigenfunctions of the Kohn-Sham equations ipeiz)'- 



2m* dz^ ~^ '^^ 



i;{z) = E^{z) (9) 



with 



1 22 7 8 -I- 4r /3 

V,^{z) = c^{z) + EMz)) + E'M^)) ^ c^{z) - ^Ry* - 0.88 / ^ Ry\ (10) 

where = {ATxn{z)a^/?))^'^/'^ ^ electrostatic energy (j){z) is related to n{z) via Eq. (j21), and 

n(z)=iV f g{e)mz)\'de. (11) 
Jo 

The detailed numerical solution of Eqs. (P)-(fTT|) will be published elsewhere. Here we 
will estimate the interface dipole term using approximate eigenf unctions of the one-electron 
Schrodinger equation with the potential in the form of a rectangular step of the height 
Ep + $s (shown in red in Fig. |21): 

mO\' = exp[-2v/l+eo-e ■ C], (12) 

where = ^s/Ep. The wavefunctions in Eq. ()12j) capture the most important feature of 
the exact solution - exponential asymptotic decay far from the interface j^]. An estimate 
for the value of the parameter Q can be obtained from the Thomas-Fermi equation inside 
n+-region (z < 0). This yields: 

Q = 0.088r;^/2ao 2. (13) 

Substituting Eqs. (fT^ and (fT^ in Eq. ((Tj) we obtain: 

3.68, 0.982 ,,,, 2.21 0.916 0.88 

*. - = ^/({„) - 7f + — + v:^s' 

where 

f( \ ^'^ 1 + 3i/ - 3yy/l + y ■ coth~^ y/l + y . . 

^^y>- 2^(1 + 3y) + {l + y)il -3y) cot-' ^ ^ 

The results of the application of Eq. (E ) are shown in Fig. IHl and compared with the exact 

numerical results by Lang and Kohn 2Q[. We see that our estimate is in good agreement 

with the exact theory for > 2.5. As we will see below this region of electron densities is 

the most relevant one for quantum computing purposes in n^/i/n~^ Si system. 

In order to apply our results to Si we have to take into account multi-valley character of 

the kinetic energy operator. If we neglect the effective mass anisotropy the main effect will 

consist in replacing parameter kp with its renormalized value (Sn'^N/ g^Y^^ , where gy is the 

number of valleys in the conduction band. The analog of Eq. ()14|1 reads: 

3.68 _ 0.982 2.21 0.916 

= ^j3-^^o - ^j^fi^o) - ^J^ + ^j^, (16J 

gv rj Ts' gv rj g^' 




FIG. 3: Workfunction of a jellium metal at different densities: blue curve - Eq. (fTl|l . circles - exact 
theory 

where the number of valleys gv=Q for Si. Here we have omitted insignificant correlation 
energy. The results of the calculations of $o by means of Eq. ()16p are shown in Fig. |3] 
We used Ry* = 34.5 meV, and = 18.3 A'24]. It can be seen from Fi g, p that for the 
purposes of quantum computing on long-lived donor states of Li in Si the electron 
concentration in n"*" electrodes must be in the range between 2.55-10^^ cm~^ {gl^^rg = 4.5) 
and 1.5-10^^ cm~'^ {gl^^rg = 2.5). The lower limit is the metal-insulator transition threshold 
for n"*" Si 3, |2^ while the upper limit corresponds to the point where the total conduction 
band offset Aq = Ep + $s between n~^- and z-layers reaches 1 Ry*. A schematic energy 
diagram corresponding to this range of concentrations is shown in Fig. El (a). It can be seen 
that the first two lowest energy levels of Li donor in Si are separated by an energy gap from 
a continuum of the many-body states in the n"'"-electrodes. Thus our qubit is protected from 
a dephasing caused by dynamic fluctuations of the electron density in the contacts. 



Qubit state readout via photo-assisted resonant tunnelling. The n^-i-n' or 
structures can be used for readout the electron states in systems considered in Refs. 
Indeed, and individual qubit in these systems is encoded on the ground and long-lived first 




Donor Concentration (cm ) 

FIG. 4: Barrier of Si n+/i junctions at different donor concentrations (blue), Eq. (|16j) . the total 
conduction band offset Aq = Ep + is also shown in red. 

excited state of a shallow donor or acceptor. According to estimates of Ref. for one of 
those systems, Li donor in Si, the hfetime of the first excited state ls,B2 (further denoted 
as |1)) can be up to 1 second if the energy separation between this state and the ground 
state ls,Ai (further denoted as |0)) is ~ 0.06 meV. The ionization energy of the ground 
state of Li in Si, /o=l Ry*=34.5 meV. The excited 2pQ-state has the ionization energy 
/p=11.5meV and the lifetime Tp ~10~^s Q]. The readout can be realized when the lower 
levels |0) and |1) lie below and the 2pQ -level is higher than (Fig.|^b)). This situation 
corresponds to doping concentration < 10^^ cm~^. In this case the lower levels of the 
qubit lie below the continuum spectrum of the contacts in equilibrium and one of them is 
populated by the electron which does not interact with electrons in conducting layers 
(Fig. El^b)). When bias voltage V ~30 meV is apphed the 2pQ -level moves down below 
the Fermi level of the left electrode and all Li levels are shifted under the action of electric 
field E = V/d^ 4KV/cm (Fig. Efc)). Radiation absorption between lower Li states |0) or 
|1) and the state 2po depends on the polarization of the infrared light. Thus, selecting the 
radiation by frequency and polarization we can excite only the state |0) or |1) depending 
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FIG. 5: Schematic band diagrams of n~^-i structure in equilibrium (a) and at a bias voltage V (b). 

on the polarization and detect the tunnelling current through the 2pQ -level of the ionized 
impurity (Fig. Efc)). The tunnelling frequency can be estimated as 

u = {EF/h)exp{-2d/3lo), 

where Iq = {h? /2m*^sY^^ is the tunnelling length. For typical values $s=8meV, 
Ep ~20meV and d=25nm we obtain uj >10^'' Hz, i.e. uTp ^ 1. It means that many 
electrons will tunnel through the excited 2pq -level before the recombination process occurs, 
the donor becomes neutral, and the current terminates. Thus a single-qubit measurement 
with the current signal much stronger than that of the single-electron measurement can be 
realized. The tunnelling current depends dramatically on Li position, namely it 
maximum when the Li atom is placed equidistantly between the electrodes 21 
effect can be used for in situ placement control of Li atoms during the growth process at high 
temperatures, when the highly-mobile Li ions are guided to their permanent final positions 
by an external electric field. 

The effect of photo-stimulated resonant tunnelling can be verified independently in the 
experiments on 5-doped Li layers placed between two n+-electrodes. These structures are 



las a sharp 
,m. This 



similar to those studied in Ref. [2^. Consider kinetics of the resonant current immediately 
after the photoionization occurs. For simplicity we assume that the tunnelling frequency is 
relatively small, u/Wdec 1, where Wdec > is a total decoherence width of the 2pQ 
level. Then the rate equations for sequential resonant tunnelling read: 

dni/dt = Ti{l -n)- r,.ni - Wni{l - uq) (17) 
drio/dt = Wni{l-no). (18) 

Here rii are the populations of the levels and 1, n = rii + uq, Ti and are inverse elastic 
tunnelling lifetimes, and W is the inverse recombination time. The first term in Eq. ()17|) 
takes into account the Coulomb blockade effect. The difference between Eqs. (fT7j) . (fT^ and 
conventional equations 30] is that the resonant current is decaying via the recombination 
channel. If Ti(^r) ^ W the population ni quickly attains quasi-equilibrium and the tunnelling 
current can be expressed as: 

I = ri{l-n)r^-^^ (19) 

where 



Therefore, the photo-stimulated resonant tunnelling current through the 5-doped Li layer 
will be characterized by a long-term decay kinetics / ~ 1/t. We note that the condition for 
the sequential tunnelling consider above is opposite to that required for the qubit readout. 
In the letter case Wt/Wdec ^ 1 and a coherent resonant tunnelling regime needs to be 
considered. 

Conclusion We proposed a device based on n~^—i—n~^ heterostructures for electronic control 
and readout of qubit encoded in two lowest states of a shallow donor of Li in Si. Qubit state 
readout is based on the polarization-sensitive photo-stimulated resonant tunnelling effect. 
A key feature of the system is that the qubit states (first two lowest energy levels of a 
donor) are separated by a finite energy gap from a continuum of the many-body states of 
the controlling electrodes. Therefore the qubit is protected from a dephasing caused by 
dynamic fluctuations of the electron density in the contacts. It is of interest to apply the 
proposed scheme for the case of qubits based on orbital states of quantum dots as well as 
on shallow acceptor states jlli] (in the later case p~^-i-p^ structures should be used). 



APPENDIX A: THOMAS-FERMI EQUATION 



According to the Thomas-Fermi theory the electrons are considered as classical particles 
obeying Fermi statistics. Then the maximal energy of an electron near point z reads: 

2 

^ma. = ^-^^<\>[z), (Al) 

where = — 6^9(2;) is the electrostatic (Hartree) energy of an electron in the field of all 
other electrons and external potential, e > is the magnitude of the elementary charge, 
is the electrostatic potential, and 'Pmax is the maximal momentum of the electron. The 
maximal energy of the electron is nothing but a chemical potential at T = 0. This chemical 
potential must be constant in equilibrium: 

Emax = const = Ep, (A2) 

where 

fe2 fc2 fc2^2 

= ^(Svr^no)^/^ = -^^(Svr^iV)^/^ = 1^. (A3) 
2m* 2m* 2m* 

Here no = is the electron density in the bulk and N is the concentration of donors. The 
maximal (Fermi) momentum can be related to the number density of the electrons n{z): 



Pr, 



KkpHz)] = h{37cMz)y/'' (A4) 



Thus, from Eqs. ()A1|) - ()A4|) the relation between electron density n(z) and electrostatic 
energy (j){z) within the Thomas-Fermi model is: 

= ^ ~ ^ (A5) 

The Poisson equation for (f) reads: 

= — (®'-^' - n) ' (^«' 

Introducing dimensionless energy u = 1 — tp/Ep-, from Eqs ()A5|l - ()A6j) we obtain the self- 
consistent, dimensionless, Thomas- Fermi equation: 

Here we introduced dimensionless length = z/Ltf, where Ltf is the Thomas- Fermi screen- 
ing length: 

/ Ef>C 1, „ o/n On /7r\l/6 1 ^ 



Here = nfi} /m*e^ is the Bohr radius. 
The equation u" = f{u) can be solved as: 



Therefore, we obtain: 



It can be shown 



—j =2 J f{u)du + Const (A9) 
22j |) that u{0) = 3/5. Therefore the solution of the Thomas- Fermi 



equation for ^ > reads: 

3 1 



u=- ,^ , ,,,,4 , (All) 



5(l + C/Co 

where Co = (lOvTS)^/^. Finally, from Eqs. ()A5|) and ()A11|) we obtain for ^ > 0: 



/ Q \ 3/2 

N 



1 + 



(A12) 



{loviEy/^LTF, 

Integrating Eq. ()A12|) from to oo we obtain the total number of electrons per unit area in 
the upper half-space (or holes in the lower half-space): 

Q^V2[^"\r.N = ^^ (A13) 

and the prefactor in the expression for 0(oo) — 0(0): 

For a multi-valley semiconductor with isotropic spherical valleys we have to replace 
Eq. (jXH) with 

Pmax = hkF[n{z)/g^] = hlSiT^niz) / gX^^ , (A15) 

where is the number of valleys. Therefore, the screening length and Q will change ac- 
cordingly, however Eq. ()A14j) will not. 

APPENDIX B: ELECTRON DENSITY 

To calculate the electron density n{z) using LDA wavefunctions we have to take into 
account the fact that we are dealing with the wavefunction of a 3D-semi-infinite metal: 

m = ^e'h%^(z) (Bl) 



The quantum number kz describes the bulk states, i.e. ipkA^) = ^"^fczl-^) ~^ Asin(A;^2; + S) 
as 2; — i> — oo. Using Eq. (IBljl we can calculate the number density: 

^(^) = TJTTTTs / / 27rA;||rfA;|,e(4 - k?, - k]) \^l^kXz)\ (B2) 

V (iTTj^ J_k^ Jo 



{kl-kl)\^PkA^)rdk, (B3) 







27r2 

Introducing dimensionless energy e = h^kl/ {2m*EF) we further obtain: 
1 /2m*Eir\^/^ /"^ 1 -e I - , ,,2 , 3^ f^l-e 



If we know Q we can normalize the wavefunctions ipe{z) in the upper half-space: 

' ' N jl deg{e) mztdz 



where 



= ^ (B6) 



APPENDIX C: ID POISSON EQUATION 



The general solution of ID Poisson equation ()A6|) reads: 

(P{z) - 0(-oo) = —— / (/ - z)p{z')dz', (CI) 

X J~oo 

where p{z) = n{z) — NQ{—z). If the system is electrically neutral, then 



0(00) - 0(-oo) = — / z'p{z')dz\ (C2) 

and 



,2 /"oo 
X 



(00) - </)(0) = / (C3) 







Substituting Eqs. ()B4j) and ()B5j) in Eq. ()(73|1 we obtain Eq. of the paper. 
APPENDIX D: MODIFIED BARDEEN TREATMENT. 

J. Bardeen, in his pioneering work of 1936 calculated a moment of a surface dipole 
layer using a model step-like potential barrier of infinite height. Her we will follow a similar 



path but assume that the height of the barrier is finite. The value of the barrier will be 
determined self-consistently. 

Let us approxiamte the exact self-consistent potential of a semiinfinite jellium metal with 
a step function: 

V,ff{z) = V{-oo) + {Ef + $.)G[L(eo) - z]. (Dl) 

Here, as previously, .^o = ^s/Ep, and we assumed that the region of the classical motion 
of electrons is extended beyond the metal surface z = 0. This assumption is necessary to 
ensure elelctroneutrality of the system. Integrating the charge densities inside and outside 
the metal and imposing the electroneutrality condition we find L{C,o)'- 



(D2) 



where 



^2(^0) 





3 




~ 16 


3 




16 


1 



+ (3-eo)cot-i(v^ 



1 + ^0 



1 + eo 

"'^''"V(l-3eo)cot-i(v^) 



(D3) 

(D4) 
(D5) 



The limiting value L{oo) = Sir/Skp was previously found by Bardeen [13] (see also Ref. 



3l|). 

As a next step we integrate the Poisson's equation and employ Eq. (|T)) to determine ^0 
self-consistently. This leads to a transcendental equation similar to Eq. (fT^ : 

0.88 



2.21 0.916 



r, + 7.8' 



where 



and 



/(eo) = 0o(eo) + 2 K(eo)-g2(eo)] 



1 + 3^0 - 3eov^r+e^coth-\v^rTe^) 



(D6) 



(D7) 



(D8) 



4(1 + eo) 

The value of obtained from Eq. ()D6jl is shown in Fig. IHl and compared with our previous 
Thomas- Fermi result. We see that the Thomas- Fermi result is in better agreement with the 
numerical calculations We believe that this is due to the fact that the total charge 



— Thomas-Fermi Treatment 

— Modified Bardeen Treatment 

O Exact Theory (Lang and Kohn, 1971) 




FIG. 6: Workfunction of a jellium metal at different densities: blue curve - Eq. (|14)) . red curve 



Eq. ()D6|) . circles - exact theory 



2^ 



within the metal is better reproduced by a self-consistent (but not oscillatory) Thomas- 
Fermi denisty rather than by an oscillatory (but not fully self-consistent) denisty obtained 
from the modified Bardeen treatment. 
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